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Abstract 

It has recently been shown that the equation of motion of a massless scalar 
field in the background of some specific p branes can be reduced to a modified 
Mathieu equation. In the following the absorption rate of the scalar by a D3 
brane in ten dimensions is calculated in terms of modified Mathieu functions 
of the first kind, using standard Mathieu coefficients. The relation of the 
latter to Dougall coefficients (used by others) is investigated. The S*-matrix 
obtained in terms of modified Mathieu functions of the first kind is easily 
evaluated if known rapidly convergent low energy expansions of these in terms 
of products of Bessel functions are used. Leading order terms, including the 
interesting logarithmic contributions, can be obtained analytically. 



-mail: manvel@phys.uni-kl.de, on leave from Yerevan Physics Institute 
-mail: muellerl@physik.uni-kl.de 

-mail:jqliang@physik.uni- kl.de, jqliang@mail.sxu.edu.cn 
-mail:ybzhang@physik.uni- kl.de, ybzhang5@yahoo.com 

1 



I. INTRODUCTION 



Recently the equations of motion of several cases of massless scalar fields propagating 
in a supergravity background describing p-brane solitons have been shown to be reducible 
to a Schrodinger-like equation with a singular potential and hence to a modified Mathieu 
equation, so that various aspects, such as absorption probabilities, become exactly calculable 
[1,2], which by AdS/CFT correspondence may yield information on correlation functions 
in a related world volume effective field theory. The singular potential appearing in the 
coefficients of the metric is in the case of the D3-brane the Coulomb potential in 6 spatial 
dimensions. In view of the fact, that very few such exactly solvable cases are known and that 
a Mathieu-type equation arises in a number of such problems as a result of the invariance 
of the wave equation under various diagonal dimensional reductions on the world volume, 
these theories are of exceptional importance and deserve to be studied in full detail. The 
two recent investigations [1,2] of the absorption of partial waves of a massless scalar field by 
D3 branes in 10 dimensions [1] and by a dyonic string in six dimensions (or a D1/D5 brane 
intersection in 10 dimensions or extremal 2-charge black hole in 5 dimensions or M2/M5 
brane intersection in 11 dimensions) [2] study the resulting modified Mathieu equation in 
terms of expansion coefficients introduced by Dougall [3] in 1916 and present very few 
details. It is not possible to follow the calculations of these papers without extensive work 
of one's own, which is made even more difficult by singularities of expansion coefficients 
that require additional attention. These studies are, in fact, complicated applications of 
modified Mathieu functions, which, in our opinion become even more complicated if instead 
of standard Mathieu coefficients, i.e. those in modern texts, the coefficients of Dougall 
are used, for the calculation of which the authors of ref. [1] developed in addition their 
own algorithm. The significance of the Mathieu equation in such contexts can also be 
seen from a different angle since the equation occurs also as the appropriately transformed 
small fluctuation equation in the study of Born-Infeld theory in the bosonic light-brane 
approximation with only an electric field E = — V0 in p = 3 dimensions and the remaining 
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components of the vector potential as massless scalar fields, reduced to only one field y in 
the simplest case. One can show that finite energy configurations of these fields independent 
of one another are not stable, but their combination with appropriate boundary condition 
[4,5] (equivalent to the Dirichlet boundary condition) is (with supersymmetry) a BPS state 
with corresponding Bogomol'nyi equation. Investigating the stability of this configuration, 
i.e. the D3-brane, with respect to transverse fluctuations of both the throat or fundamental 
string and the brane, one again arrives at an equation with the singular potential 1/r 4 [4,6] 
which can be converted into a modified Mathieu equation [7,8]. Thus in each of these cases 
a Schrodinger-like equation is obtained with the singular potential 1/r 4 . Such potentials 
have been the subject of investigation 30 years ago [9] and were motivated by the lack of 
understanding of weak interactions at that time. Thus the potential 1/r 4 and the associated 
scattering problem had also been investigated, and various forms of the S-matrix had been 
given [7,8,10,11] in terms of modified Mathieu functions or related functions for which - at 
the latest since the publication of refs. [12] and [13] - widely used definitions and notations 
exist. 

In view of the scarcity of fully solvable examples of theories on a supergravity background 
we consider it worthwhile to reexamine the case of the propagation of a massless scalar field 
in the presence of a 3-brane by using modified Mathieu functions with standard Mathieu 
coefficients and the S'-matrix evaluated in terms of these. In our opinion these calculations 
are more transparent than those using Dougall coefficients and are easier to follow with 
reference to modern literature on the subject. In view of the complexity of the calculations, 
due also to the fact that later iterations contribute to earlier lower order terms, we present 
these in some detail. Our presentation below should therefore also enable others to follow 
the reasoning, and this particularly since leading order terms can be understood without 
resorting to numerical methods. 

In the following we first formulate the semiclassical gravity problem and reduce it to the 
modified Mathieu equation. We do not rederive the S'-matrix, but recapitulate in Appendix 
A the main steps in the derivation, and in particular some steps that have not been written 
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out explicitly in ref. [8], this being the prime reference on which our considerations are 
based. We then consider briefly the gauge field theory approach in a simplified Born-Infeld 
version to demonstrate how this also leads to the Mathieu equation. Following this we 
consider the Floquet exponent associated with Mathieu functions and show how this has to 
be calculated in singular cases (such as the the cases to be considered here and in the S- 
wave case already in the dominant approximation). The calculation of coefficients of series 
expansions of modified Mathieu functions is then considered and the Dougall coefficients 
used in refs. [1,2] are compared with ordinary, i.e. standard, Mathieu coefficients as in ref. 
[12]. We calculate examples in singular and asymptotic cases (the latter being those that 
permit one to ignore the singularities of early coefficients). Higher order contributions are 
obtained with Mathematica. We show that Dougall coefficients are more difficult to obtain 
than ordinary coefficients - an observation that may explain why Dougall did not evaluate 
any of his own coefficients in his work of 1916. We then evaluate the relevant quantities 
appearing in the S'-matrix and hence the absorption probabilities and cross sections. Where 
comparable, our results can be seen to agree with those of ref. [1]. The treatment presented 
below makes full use of the well established theory of the Mathieu equation and can therefore 
point the way to explore other aspects, such as application to double-centered D3 branes 
and to higher energies which have been discussed recently [14]. 

II. THE SCALAR FIELD IN THE D3-BRANE METRIC 

The supergravity background for an extremal Dp-brane in the 10-dimensional type IIB 
theory is [15-17] 



where (r being the radial coordinate in the SO (5) symmetric space orthogonal to the branes) 




(1) 



v 




(8-P) 
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i=i 
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and the harmonic function H is given by 



H = l + 



(3) 



For p — 3, the case of interest here, i.e. the D3 brane coupled to the 4-form RR-potential 
[17,18], R with R 4 = 4ng s Na' 2 (g s the string coupling and iV the number of D3 branes) is 
the radius of S 5 and AdS 5 in the socalled decoupling limit in which one obtains a duality 
between M — 4 U (N) supersymmetric Yang-Mills theory in 4 dimensions and string theory 
in the near horizon AdS 5 x S 5 background [19,20]. As pointed out in ref. [1], for a comparison 
of considerations in terms of supergravity and those in terms of D-branes, one is interested 
in the domain of small uR, where u is the energy of the field incident on the brane. 
For a massless scalar fluctuation field around the dilaton field <3> given by [15] 

(which is constant for p = 3) in the background of this metric, the equation of motion is 



-pd^y/gg^d^ = 



(4) 



After separation of the S 5 harmonics Y(Qj), in particular the Gegenbauer polynomial 
C;(cos#), where x = rcosO, and a factor e tult the radial wave function ipi(r) = y(r)/r% 
of the l-th partial wave of energy u of the scalar field is found to satisfy 



r 5 dr \ dr 



1(1+4:) 



u 2 R A 



d 2 



+ oo z + 



2 i? 4 (/ + !)(/ + 



y = 



(5) 



^dr 2 r'* r z 

We see that for R 4 ^ this is the equation of an attractive singular potential with coupling 
constant g\ = u 2 R 4 . For u 2 > an incident wave allows both transmitted and reflected 
waves, and from the ratio of coefficients one can determine the ^-matrix. It is convenient 
to make the substitutions 



y = r 1 / 2 0W, r = 7 e 2 , *y = g /h, K 1 = u 9o = u 2 R 2 , A = (I + 2) 2 , 



(6) 
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+ 



2h 2 cosh 2z - A 



= (7) 



which convert the range of r from to oo to that of z from — oo to +oo. The equation 
thereby becomes the modified Mathieu equation 

dz 2 

In view of the principal interest in the relation of our semiclassical gravity consideration 
with the superconformal limit of the dual theory in the near-horizon domain, we are here 
interested in waves of low energy, i.e. of small u, and so in solutions of the modified 
Mathieu equation around h 2 = 0. The modified Mathieu equation allows series expansions 
of this type in terms of exponential, hyperbolic and cylindrical functions and (surprisingly) 
in each of the cases with the same coefficients C2 r {h 2 ) where v is the Floquet exponent and 
the subscript r a positive or negative integer or zero (not to be confused with the radial 
coordinate). The solutions in terms of exponentials are written Me u (z, h 2 ), those in terms of 
hyperbolic functions cosh and sinh Mc u and Ms u . The solutions of the i-th kind are those 
in terms of cylindrical functions and are written Mu\z,h 2 ) where % = 1,2,3,4 correspond 
respectively to expansions in terms of Bessel, Neumann or Hanked 1,2 ) functions. The series 
of M$\z, h 2 ) converge uniformly only for | coshz| > 1, whereas the series of Me u (z, h 2 ) is 
uniformly convergent for all finite complex values of z as shown in ref. [12]. Since r = 
corresponds to z — — oo, the S-matrix is obtained by continuing the solution at z — — oo to 
z = oo. This means that a solution has to be continued, via matching to Me± v (across 
the domain \z\ < 1), to a linear combination of MjP and at +oo. A few main steps of 
this calculation are given in Appendix A. The expression for the S'-matrix finally obtained 
is 

where 

Ari"(0,A 2 ) 

This S'-matrix describes the scattering of an incident wave (component of the scalar field) 
of energy uj off the spherically symmetric potential. One could visualise this scattering as a 
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spacetime curvature effect or - with black hole event horizon zero - as that of a potential 
barrier surrounding the horizon. With the horizon shrunk to zero at the origin (implying 
in the field theory a relation between mass and charge reminiscent of the Bogomol'nyi 
equation), this extremal case corresponds to that of a BPS state. 

The absorptivity is A = 1 — SS*. The absorption cross section differs from this by a 
multiplicative factor in front. The absorption cross section a l abs of the l-th partial wave in 
n spatial dimensions has been derived in ref. [21] and is given by 

i . 2^-V*/ 2 - 1 . /o .,. ; , /o /Z + n -3V ,_ |2 



<J abs 



D 1 



_(n/2-2)!(Z + n/2-l)^ + \ 'j(l-|5| 2 ) (10) 



For n = 6 as in our case this /-wave (here semiclassical) absorption cross section (or socalled 
greybody factor) is given by 

8tt 2 



a abs — 



3u 5 



(/ + 1)(Z + 2) 2 (1-|S| 2 ) (11) 



III. THE D3 BRANE IN BORN INFELD THEORY 

To supplement the previous section, we consider briefly the simplest version of super- 
symmetric Born-Infeld electrodynamics for the 3-brane. Our main intention is to recall 
that the equation of small fluctuations about the _D3-brane is again a modified Mathieu 
equation as obtained above. In the simplest such model reduced to the static case we write 
the Lagrangian 



d p xC, £ = 1 



1 - {d^f + («%) 2 + 



- £ p e<M(r) (12) 



Here E i = F 0i = —di(p,i = 1, ...,p and y(xi) originates from one of the gauge field components 
A a for a = p + 1, ...,d — l,d= dimension, which represent the transverse displacements of 
the brane (of which we consider only one, e.g. A 9 ). The source term of the electric field 
with charge e and S 3 = 4n hints at spherical symmetry. Considering only this case here and 
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hence that of S'-waves, we obtain two Euler-Lagrange equations which we can write 

dr {^mi)) = °' ^Wv) = c 

where c is a constant of integration. Explicitly, 

0' e —y' c 



[i - (0O 2 + " rP ~ v [i - {<t>T + (y') 2 ] 1 rP_1 

so that 



(13) 



£ _ e = 1 
y' c a 



(14) 



Then 



(^') 2 r 2(p-l) + e 2(l _ fl 2) ' ^') 2 r 2(p-l) + e 2(! _ fl 2) ( 15 ) 

The p-brane and anti-p-branes are now given by 

+ poo y 

y(r) =(— ) ae / rfr — (16) 

A ^/ r 2(p-l) _ r 2(P-l) 

where r 2( - p ^ = e 2 (a 2 — 1) > 1. In view of the proportionality (14) the Lagrangian can be 
written 

£ = 1 - y/l-(l-a?)(di<l>y - S p e05(r) (17) 
The contribution to the energy not including the source term is for p = 3 



= / Sx{ - l = 
J \ v/1 - (1 - a 2 



£„, = / (fx { = = - 1 

-(l-a 2 )(W 

Only for a charge e which is kept fixed under a scale transformation is the energy minimal 
in the limit a 2 — > 1. This is the limit of the Bogomol'nyi bound and hence for this value 
of a 2 the Born-Infeld configuration, i.e. the Dp-brane or string is classically stable, i.e. a 
nontopological BPS state. The reason is that for a = 1 we have 0' = y' which in the original 
context with y = Ag implies F 0r = d r A 9 . This again has implications for the supersymmetry 
variation of the gaugino Xi the susy partner of the gauge field, which is 

5 X = T^F^e 



(//, v being the original 10 dimensional indices and Y^ v the appropriate combination of 10 
dimensional Dirac matrices). In the case of only the electric field and the one excitation 
under consideration the gaugino variation is 

5 X = -(T 0r d r A + T 9r d r A 9 )e 

a ^ -(T 0r + T 9r )d r A e (18) 

Thus, since 1 + r°T 9r is a projection operator it is precisely for a = 1 that the variation 
5x can be zero for some nonzero e, thus preserving a fraction (half) of the number of 
supersymmetries. 

The energy E ns is infinite, but integrating from r = 5 to infinity so that y = 
ae/S, dE ns /dy, the energy per unit length of the string is finite, i.e. 

dE ns 1 2 4vre 

— - — = -(1 — a ) = const. 

As shown in [4], unless a = 1 supersymmetry is completely broken (i.e. the supersymmetry 
variation of the gaugino would not be preserved). In this limit the throat radius r becomes 
smaller and smaller, and the brane pair moves further and further apart. If one considers 
small fluctuations £ orthogonal to the brane and the string, one obtains the small fluctuation 
equation [4] 



1 + 



r 2(p-l) 



e = o (19) 



where f2 2 > for stability. The radial part of these equations is with £ = r (P 2 and 
angular momentum I 



d 2 j) 
dr 2 



+ 



1 /,„ , „ (P-1)(P-3)1 , o2 A , e 2 (p-2) 



Z(Z + p-2) - ^ ^ ^ ^ + fi 2 1 + 



\2 



ip = 



r 2 ] v ^ 7 4 ( 1 y r 2(p-i) 

Thus for p = 3 and x = fir, « = ef2 2 and for S-waves since the string cannot depend on the 
angular variables of the worldvolume 

(£ + 1 + ^ = (20) 



\ ofe 2 ar 



This equation is an S-wave radial Schrodinger equation for an attractive singular potential oc 
x~ 4 but depends only on the single coupling parameter k with constant positive Schrodinger 
energy. 



IV. CALCULATION OF THE FLOQUET EXPONENT IN SINGULAR AND 

NONSINGULAR CASES 

The Floquet exponent v enters the discussion of the Mathieu equation in view of the 
Bloch wave property of the modified Mathieu function Me v (z, h) 

Me u (z + in, h) = e iU7T Me u (z, h) (21) 

The Floquet exponent can be introduced in several ways. In ref. [12] (p. 107) v is introduced 
by the relation 

cosnu = yi(n, A, h 2 ) (22) 

where yi(x) is a fundamental solution of the (periodic) Mathieu equation satisfying the 
boundary conditions yi(0) = l,y'j(0) = 1 and A is the eigenvalue which, of course, is not 
necessarily an integer. With a perturbation theory ansatz for yi(ir, A, h 2 ) around h 2 = the 
following expansion is then shown to result (cf. ref. [12], p. 124): 

r- 4 TTSUITiVA 
COS IIV = COS7TV A + h — = h 

4VA(A- 1) 



+ h b 



15A 2 -35A + 8 r- 7r 2 cos7r-\/A 

-n sin 7tv A — 



+ 0(h 12 ) (23) 



.64(A- l) 3 (A-4)Av / A 32A(A-1) 2 _ 
Alternatively one can apply directly perturbation theory to a trivial periodic solution of the 
limit h 2 — > as shown in ref. [8]. In this case the following expansion is obtained 



h 4 (5z/ 2 + 7)h 8 



A = IT + 777— n 7T + 



+ ^#^^ + 0(O (24) 



2(z/ 2 -l) 32(z/ 2 - l) 3 (z/ 2 - 4) 
9z/ 4 + 58z/ 2 + 29)h 12 
64(z/ 2 - l) 5 (z/ 2 -4)(z/ 2 - 9) 

This series can be reversed to yield u, i.e. 
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, // ; (13A-25)/i 8 



z/ = A - 



2(A-1) 32(A- l) 3 (A-4) 
(45A 3 - 455A 2 + 1291A - 1169) h 
64(A- l) 5 (A-4) 2 (A-9) 
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+ 0(h 16 ) (25) 



and so 



„ = Vx + * ^ - (8 - 35A + 15A2) ^ + 0(0 (26) 
4(1-A)VA 64(A-4)(A-1) 3 Av / A 

One can easily check the agreement with the above expression for cos nu by evaluating cos nu 

with v of this expansion. 

An obvious feature of all of these expansions is that they are singular for integral values 

of A. This behaviour is wellknown. It means that these expansions are in these cases really 

asymptotic expansions for large values of v or A and can be used in such cases. However, 

for other values the expansions can also be convergent for sufficiently small values of h 2 as 

is shown in ref. [12]. For values of v close to an integer or \f\ close to an integer one has to 

expand around these as is also mentioned in ref. [12] (pp. 124-125). We demonstrate this in 

the case of v almost equal to 2. Thus we set 

v = 2 + 5, \ / \ = 2 + e 

so that 

7l 2 5 2 

cos ixv = cos 7r(2 + 5) = cos 2n. cos n5 = 1 h • • • (27) 

and consider the limit e — > 0. Expanding the cosine and sine expressions appearing in eq.(23) 
about A = 4 we have (A — 4 fx 4e) 



and 



nr n . w . nr N 7T (A — 4) 2 7T 2 

cos y \n = cos27r + A — 4 )( — sin v A7rh =4 . — = + • • • = 1 

2V\ 8 



sin V\tt = sin 2tt + (A — 4) (cos V\tt) a=4 . — = + • • • = — - — + 

2^ 2^ 



+ 



Substitution into eq.(23) and considering the approach A — > 4 gives 
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COS 7TU = 1 (A - 4) 2 + • • ' 



4v / A(A - 1).2\/A 
„/ (15A 2 -35A + 8)tt 2 (A-4) tt 2 .1 \ 

+ V64(A-l) 3 (A-4)Av / A.2v / A ~ 32A(A - 1) 2 J + 
1 (A-4) 2 vr 2 \ h^ 2 ( I \ A y 

+ 100 ,ty iv* f( 11A ' - 31A + 8) + • • + 0(/i 12 ) (28) 



128A 2 (A - l) 3 

Hence (observe the cancellation of factors (A — 4) in the term of 0(h 8 ))'m the limit e — > 

, /i 8 7r 2 (ll x 16-31 x 4 + 8) , 5vr 2 /i 8 

«)8 7ri/ = l + ^ - + •■• = ! + -293T + --- ( 29 ) 

Comparing the expansions (27) and (29) we obtain 

We see that although the coefficients of higher order terms of eq.(23) contain factors (A — 4) 
in the denominators and so suggest divergences, the trigonometric factors sin7r\/A in the 
numerator always cancel these out and thus yield a regular expansion for v which is even 
convergent within a certain domain around h 2 = 0. One thus obtains the expansion 

"= 2 -¥w 2 ) 4+ i^^ + S (k/2)a - (30) 

This expansion has also been given in ref. [1]. Expansions around other integral values of 
are obtained similarly. Of course, the higher the value of this integer, the more terms at 
the beginning of the series are identical with those obtainable from the perturbation series 
(26) above. Thus in the case of \/X = 3, 4, we obtain from the first two terms of expansion 
(26) 



and similarly 
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in agreement also with results of ref. [1] up to the given order, i.e. 
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v/A = 3 : u = 3- -(h/2f + —(h/2) 8 + — (h/2) 12 + ■■-, 
6 V ' ' 4320 V 1 ' 155520CT 1 ' ' 

V\ = 4 : v = 4 - -(V2) 4 - ^^(V2) 8 + 3 ° 5843 (h/2) 12 ■ ■ • (31) 
15 l 1 ' 27000 1 7 ; 680400000 1 1 ' y ' 



V. RELATION BETWEEN DOUGALL AND STANDARD MATHIEU 

COEFFICIENTS 

One may wonder how the Mathieu function coefficients of Dougall [3] which are used in 
refs. [1,2], are related to those in modern standard literature such as ref. [12]. We therefore 
demonstrate their precise connection here. It is crucial thereby to distinguish between 
nonsingular or asymptotic cases and singular cases, as we shall see. We begin with the 
nonsingular case and calculate a coefficient given in ref. [12] (up to the first nonleading 
contribution) by starting from Dougall's definition of his coefficients. We shall see that the 
coefficients given in ref. [12] obtained from simple continued fraction solution of the basic 
recurrence relation of the coefficients are in this case not only easier to derive but have also 
a simpler form than the coefficients of Dougall. 

The modified Mathieu function in terms of exponentials is defined in ref. [12] as the 
following sum 

oo 

Me u (z,h 2 ) := cUh 2 )e {u+2r)z (32) 

r=— oo 

where v ^ ±1, ±2, • • •. In ref. [12] (p. 131) the following relation of general validity is given 
and used 

Me_„(z, h) = Me u (-z, h) (33) 

This relation implies that 
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We shall see explicitly that this relation holds also in what we call the singular case below. 
Dougall [3] defines in his work the solution corresponding to Me u (z, h) as 



J(u, z) 



We therefore expect the equivalences 

- Me (z h 2 ) ( ir 0(n + i//2) - (36) 

We now verify the latter of these relations for the case n — 1 in the nonsingular case (i.e. 
for v 7^ integer+0(/i 2 )), i.e. we show that 

0(^/2 + 1) <? 2 {h 2 ) h 2 {u 2 + 4 v + 7)h 6 

<j>(v/2) ~* c£(/i 2 ) 4(i/ + 1) 128(z/+l) 3 (z/ + 2)(z/-l) + 1 ' 

where the expression on the right is given in ref. [12] (p. 121). We also show thereby that in 
leading order for small h 2 

0W2) -(n + ^)!n! (1 + ° (M) (38) 

in agreement with ref. [12] (p. 121). The demonstration of agreement requires eqs. (24), (25) 
(cf. also [12],p. 119), i.e. 

, o 2 2 & 4 (5z/ 2 + 7)/i 8 

S ^ + 2 ' « =" + 27^1) + 32(, 2 -l)3(, 2 -4) + --- (39) 

or 

±s = u + —— + • • • (40) 

4z/(z/ 2 - 1) 

This general relation is an asymptotic expansion in v (i.e. for v large), and can be obtained 
perturbatively [8]. It is crucial, of course, to deal separately with values of v close to a 
singular value like v — 2 (see below). 

Dougall defines his coefficients by an expansion, of which the leading and next-to- 
leading contributions are [3] 
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(f>{n + u/2) 



E 

Pi=0 



(h/2) 



2n+u 



(n + u/2 + s/2)\(n + u/2 - s/2)\ 

my 



1 - 



[n + v/2 + s/2 + l+ Pl ) 



1 



(n + u/2 + s/2 + 2 + Pl )(n + u/2 - s/2 + 1 + Pl )(n + u/2 -s/2 + 2 + Pl ) 



+ 



Taking into account only the leading contribution, we have 

0(n + u/2) (h/2) 2n (u/2 + s/2))\{u/2 - s/2)! 



4>{v/2) 



(n + u/2 + s/2)\{n + u/2 - s/2)! 



Using 



(-z)\ 



7T 



(z — 1)! sin7T2; 
and the approximation s u (cf. eq.(40)) we obtain 

0(n + u/2) _ {h/2) 2n u 



0(h 4 ) 



(41) 



(42) 



(43) 



(44) 



<P(u/2) (n + u)\n\ 

in agreement with c y 2n {h 2 ) / '^(h 2 ) of ref. [12]. We see therefore that the expansion (40) plays 
an important role in establishing the connection between the coefficients 0(r + u/2) and <f 2r 
in the nonsingular case. 

We now consider the next-to-leading contribution in eq.(41). Setting 

1 



W = E 



Pi=0 



(q + s/2 + 1 + Pl )(q + s/2 + 2 + Pl )(q - s/2 + 1 + Pl )(q -s/2 + 2 + Pl ) 



(45) 



we can make the partial fraction separation 

1 ( 1 



= E 

pi=0 



1 



1 



s(s - 1) \ Pl + g + 1 + s/2 pi + q + 2 - s/2 
1 1 



(46) 

s(s + 1) \ Pl + q + 2 + s/2 Pl + q + l-s/2j\ 
The sums over individual terms are divergent. Thus the expression depends crucially on 
taking differences. In order to deal with these we use the formula [22] 



n— 1 

~T1 



k=0 



k + y 



i){n + y) -ip(y) 



(47) 
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where ip(y) is the derivative of the log of the gamma function T(y). Thus 

n-l 



E 

Pi=0 

and 

n-l 

E 

Pi=0 

so that 



1 



Pi + g + 1 + s/2 pi + g + 2 - s/2 



^(g + 2-s/2) - ^(g + l + s/2) 



Pi + q + 2 + s/2 pi + g + 1 - s/2 



V>(g+l-s/2)-^(g + 2 + s/2) 



Al 1 > = 



s(s-l) 
1 



+ 2 - s/2) - i){q + 1 + s/2) 



^j{q + 1 - s/2) - ^j{q + 2 + s/2) 



s(s + 1) 

We now use eq.(40) in order to reexpress s interms of v. Then 



(48) 



1 



s(s-l) v(v-\) 

1 1 



s(s + l) + 



1 - 
1 - 



{2v- l)h 4 



4u 2 (u- l) 2 (u + l) 
{2u+l)h A 



+ ■ 
+ ■ 



(49) 



4v 2 (v- l)(v+ l) 2 

Setting q = n + u/2 and dealing similarly with the arguments of the functions ip, we obtain 
in lowest order of h 



4 (1) = 1 

q=n+u/2 _ ^ 



Again we consider a difference, i.e. 



ip(n + 2) -tp(v + n + 1) 
ip(n + l) -ip(u + n + 2) 



(50) 



A„ := A. 



(i) 



We now use the formula 



n+u/2 n-l+v/2 



^(n+1) = -C + 1 + - + I + ... + - 

2 6 n 



(51) 



where C is the Euler constant. Then 



il>(n + 2) -tp(n+ 1) 



1 



n+1 
16 



, ip(n + 1) - ip{n) = - 



and A n becomes 



1/(1/ -1) 
1 



n + 
1 



y - + n + 1) + + n] 



n 



- i>{v + n + 2) + ^(i/ + ra + 1) 



i/(^ + l) 

We now require yet another formula of the ip function, i.e. 



(52) 



#*) = -c + E(^t-— 1 

v 7 ^ \n+ 1 x + n/ 



(53) 



With this we obtain 



t])(v + n) — ip(v + n + 1) = 



-, ^(z/ + n + l) -^(z/ + n + 2) 



(where in each case the dummy summation index of the second sum, i.e. i, was renamed 
i — 1). We therefore obtain in the dominant approximation 



For n = 1 this implies 



1 



n + 1 v + n 



n v + n + 1 



(54) 



Al _ 4(1) _ 

^1 _ A q=l+v/2 A q=v/2 



1 



(55) 



2(i/+l)(i/ + 2) 

in the dominant approximation. This difference will now have to be substituted into the 
Dougall coefficient 



<t>{ 



, /2 + 1} (|) 2 W2 + */2)!(i//2 - */2)l{l - (|)M« +1 + ■ ■ ■} 



0(Z//2) (*/2 + 1 + s/2)\(v/2 + 1 - S /2)l|l - {\YA% + ■ ■ j 

,x 2 {l-(l) 4 (^ +1 -^)+---} 
7 (z//2 + l + s/2)(z//2 + l-s/2) 



Inserting eq.(55) and for s the series of eq.(40) we obtain 
0(z//2 + l) h 2 h 6 {u 2 + Au + 7) 



(56) 



+ 



0(i//2) 4(z/+l) 128(i/ - l) 3 (^ + 2) 



+ 0(/i 



10 \ 



(57) 
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in agreement with ref. [12] (p. 121). To obtain the Dougall coefficients in this form is thus 
seen to be rather complicated. This may explain why Dougall himself does not evaluate any 
of his coefficients explicitly. 

In the remainder of this section we calculate with the method of Dougall the important 
coefficients <fi(±v/2) for the S-w&ve case, i.e. s — 2, and demonstrate the agreement with 
results of ref. [1]. The results will also be needed in the next section in establishing ratios 
corresponding to the ratio of eq.(57). From eq.(48) we obtain for the leading term in the 
limit h 2 — * for q = v/2 and s = 2 (i.e. v 2, a socalled singular case) 



A 



(i) 

v/2 



i) (v/2 + 1) - i>{yj2 + 2) 



i>{yj2) - tl)(y/2 + 3) 



h 2 ^0 



V(2)-V(3) 



Using eq.(51) we obtain 



A {1) - — 



(58) 



The authors of ref. [1] developed another algorithm in which (cf. their Appendix A) 



A (i) = A (i) _ 



(59) 



and for s = 2 (their r = 1) 



3 + 2z ip(z) 



-4>(z + 2) 



3z(2 + z) 



(60) 



Here 



^z)-^(z + 2) = ^- 



hiY{z) -\nT{z + 2) 



2z + l 
'z(z + 1) 



(61) 



so that 



A^ 



3z(z + l)(z + 2) 



(62) 



and for v 2 one obtains again 



1 



V2 lg 



as above. It follows that for this case with $tv > (cf. eq.(41), 



0(^/2) 



(h/2) 




(V2) 4 



(is/2 + s/2)\(v/2- s/2)\ 



\{h/2f l + 0(h 4 ) 



(63) 



Fortunately for this case the expansion (57) does not seem to possess terms which diverge 
for v close to a positive integer^ 1, and evidently even if it did - since in our cases v = 
an integer ^ 1 + 0(h 4 ) - would not affect the leading term of eq.(59). This is radically 
different when v is close to a negative integer such as —2 in that case. We see from (53) 
that if x = — 2 + 0(h 4 ), tjj oc l/h 4 , and hence oc l/h 4 , and so this term will contribute 
to the leading factor in the coefficient of eq.(41). 

We now consider the coefficient 0(— v/2), 9fa/ > 0, for which 



and calculate this first with the method of Dougall. The following steps given explicitly 
demonstrate clearly how singularities in the limit ti 2 — > arise and how they have to be 
handled. Thus with eq.(41): 



s = 2, v = 2- 




(64) 
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<j>{-v/2) 



{h/2)- 



-u/2 + s/2)\(-u/2-s/2)\ 

(V2) 4 



E 

Pi=0 



-v/2 + s/2 + 1 + pi)(-i//2 + s/2 + 2 + Pl ) 



-z//2 - s/2 + 1 + pi)(-u/2 - s/2 + 2 + Pl ) + 

(V2)- 2 



0!(-2 + ^(/i/2) 4 )! 

1 



1 - (ty2)< 



(1)(2)(-1 + ^(V2) 4 )(^#(V2) 4 ) 



+ 



1 



(2)(3)(^(V2) 4 )(1 + ^(V2) 4 



pi =2 



(ty2) 



-2 



(-2 + ^(V2) 4 )! 



(V2)" 



(-2 + ^(/i/2) 4 )! I ^ 



1 - (V2) 4 
2 



+ 



-2)(^(V2) 4 ) 6(^(V2) 4 ) 



+ 



1 + 



eq 



(43) 1 



7T 



(V2)" 2 1 + t 



"(V2) 2 



z^5 
2 + i\/5 
6 



;i)(-l)vr^(V2) 4 



(65) 



The steps above clearly show how the singular terms in the next -to-leading contribution 
contribute to the dominant order. 

The result (65) will now be shown to agree with the calculations of the method of ref. 
[1]. For this purpose we set in eq.(62) z = —v/2 and replace v by the expression in eq.(64). 
Then 



L (i) 

l -u/2 



-8 



3.2(«v^/3).2.(/i/2) 4 iy/h{h/2Y 



and 



^i-(V2)M«=i±^ 



From eq.(41)we obtain therefore 

4>{-v/2) 



my 



0!(-2 + ^(/i/2) 4 )! 



.v 



(66) 
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Using for the factorial again eq.(43), one obtains 

7T 



0(-z//2) = -^^ 6 W ; J l + _ - I (67) 



Thus finally 



ft-"/*) = - 2 ~^{h/2f (68) 



in agreement with eq.(65). 



VI. CALCULATION OF STANDARD MATHIEU COEFFICIENTS 

Our next objective is to compute a Dougall coefficient (i.e. a ratio of two quantities 0) 
and to compare it with a standard Mathieu coefficient in the nontrivial singular case. As 
a suitable example we choose the coefficient <fi{—v/2 + 1) which when divided by <j>(— v/2) 
(calculated above) ought to agree with the Mathieu coefficient —c 2 l/ /c Q u according to eq.(36), 
i.e. we wish to demonstrate that 

<f>(-v/2 + 1)M-*V2) = -cjVcb " (69) 

We begin with the calculation of the Dougall coefficient. Using eq.(41), we have 

(h/2f- v 



-v/2 + 1) = 



(l-z//2 + s/2)!(l-z//2-s/2)! 

(V2) 4 



()V .. z//2 + s/2 + l+p 1 )(l-z//2 + s/2 + 2+p 1 ) 

1 

(1 - v/2 -s/2 + l+ Pl )(l - v/2 -s/2 + 2+ Pl ) + " 

eq.(&A)forv 1 



+ 



6 

1 



1 - (h/2) 4 



(l)(-l + ¥(V2) 4 )!l ^' " L(2)(3)(^(V2) 4 )(1) 

oo 



(3)(4)(1)(2) 



pi =2 



(0)!siii7r(l - f x _ 1 



7T [ iyfh, 

(h/2)' (70) 



6 
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With eq.(68) we obtain therefore 

0HV2 + 1) 

<t>{-v/2) 

Similarly one obtains 



(h/2f 



1-iVE 
2 + i^5 



l + 0(h 4 ) 



(71) 



0OV2 + 1) = -(V2) 4 



l + 0(h 4 ) 



+ V - h h/2) 2 

<Ku/2) ~3 W2} 



l + 0(h 4 ) 



(72) 



Our next step is to derive the corresponding expression from the continued fraction 
relation of the recurrence relation of the standard Mathieu coefficients. This recurrence 
relation is given by (cf. ref. [12], p. 117) 



'2r 



1 



(73) 



2 -- 2 h- 2 [s 2 - (z/ + 2r) 2 ] - 



h~ 2 [s 2 - {u + 2r + 2) 2 ] - 



1 



h~ 2 [s 2 - (i/ + 2r + 4) 2 ] • •• 
Here we set r = 1 and replace v by —v. Then we again use (64) and a) set s — 2, and b) 
replace v by the expansion given in eq.(64). One then has to go as far as the terms explicitly 
written out in the following continued fraction only to obtain the dominant contribution: 

1 



"2 
— v 



(74) 



/TV -(-" + 2)1 - 



^-2 [s 2_ ( _ i/ + 4)2] _ 



1 



h~ 2 [ s 2 -(-u + ey 



Making the substitutions we obtain 



h- 2 [A] 



(75) 



: (V2) 4 ] - 



/r 



■12] 



This can be seen to reduce to 



% = (h/2f- 



(76) 



' o iv^5 + 2 

which agrees with the negative of the above Dougall coefficient as expected on the basis of 
eq.(35). One can see that the calculation here is simpler than that of both (68) and (70). 
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The reciprocal of the continued fraction relation (73) is (cf. ref. [12], p. 117) 



c 



2r-2 



"Ir 



h- 2 [s 2 - (u + 2r-2Y 



h- 2 [s 2 - (z/ + 2r-4) z 



h- 2 [s 2 - (u + 2r-QY 



(77) 



Replacing here v by — v we have 



-2r-2 



-2r 



h' 2 [s 2 -(-u + 2r-2) 2 ] - 



h- 2 [s 2 - (-u + 2r - A) 2 ] - 



h- 2 [s 2 - (-u + 2r-QY 



Here we put r = and again make the replacements of eq.(64). Then 



^ = 12 



l + 0(h A ) 



For r = 2 in eq.(78) we obtain 



12 



<% v h 2 {l-iV5) 



so that with eq. (76) 



(78) 



(79) 



(80) 



— V —V —V o 

Cq c 2 c o 2 + ^V5 



In a similar way we obtain 



c'l h 2 



--2 



2 5 



l + 0(/i 4 ) 



C ^ C ^ C 2 ^ 



Cri C 9 Cn 3.2 



l + 0(/i 4 ) 



(81) 



(82) 



-0 "--2 <-0 

Summarising we have as leading contributions of standard Mathieu coefficients in the 
singular case / = or s = 2: 
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= (h/2) 



-2_ 

v 



—v 
4_ 

v 



— v 
4 
v 



— v 



iy/E-1 



{h/2? 



l + 0(h 4 ) 



2 + i^/h 
(h/2? 



l + 0(h A ) 
1 + 0(h 4 



2 3 .3 
(h/2) 2 



1 + 0(h 
1 



2 + i\/5 . 



l + 0(/i 4 ) 



{h/2f 



l + 0(h 4 ) 



(83) 



c " 2 3 .3 2 .5 

In Appendix B we give several more terms of these expansions calculated with Mathematica. 



VII. EVALUATION OF THE QUANTITY R IN THE SINGULAR CASE 

Having determined the standard Mathieu coefficients in the singular S-wave case, we 
can proceed to evaluate the quantity R entering the S-matrix. R was defined in ref. [8] (see 
also Appendix A) as 

R = a u /a„ u , a u (h) = Me u (0, h)/M^(0, h) (84) 

The function Me u (z, h) was defined previously. The functions M$ (z, h) , for % = 1,2, 3, 4, are 
corresponding expansions of the modified Mathieu function in terms of cylindrical functions 
J u (z),Y u (z), H^\z), Hl 2 \z) respectively. In particular we have the expansion (cf. ref. [12], 
p. 178) 

oo 

Me v (0,h)Mj > 1 \z,h)= c u 2r (h 2 )J„+2r(2hcoshz) (85) 

r=— oo 

As shown in ref. [12] (p. 180), a much better expansion to use in practice for M^\z,h) in 
view of its rapid convergence, is 

+oo 

4^ 2 )mS(z, h)=Y, (-l) l 4r(h 2 )J^ r (he- z )J ±l/+l+r (he z ) (86) 

l=— oo 
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so that in particular 



+00 



(87) 



i=— 00 



This formula is amazing. It implies that one and the same quantity M±l(0,h) can be 
obtained from many different expansions (and so different Bessel functions) by allocating 
different values to r, i.e. e.g. r = and 2. An analogous observation has also been made by 
Dougall [3]. 

We begin with the evaluation of Me v (0 : h), i.e. 

cUh 2 ) 



With the help of the standard Mathieu coefficients evaluated previously we obtain 



(88) 



Me„(0, h) 



Co(0) 



I - 



2 + iV5, 



Co(0) 



iV5 



2 + 



Me_ u (0,h) 



(89) 



h 2 ^0 



The last equality follows also from the general property Me- V (z, h) = Me u (—z, h) (cf. ref. 
[12], p. 131). More terms can be calculated with Mathematica. Thus 

m rn h) 1 i c * {h2) i c ^ h2) i c ^ h2) i c ^ {h2) i c ^ {h2) i c ^ {h2) i c - s(/l2) i c - wih2) 



c" (/i 2 ) 

We have to take into account also ~ 8 ,, 2 ^ 
and h 6 , i.e. 



and ~l°( h2 \ because these contribute to orders h 

c (n ) 



VE) fh\ A 202^ + 35^ fh 



(h 2 ) 2 3 3{t + V 2 / 2 4 3 3 5(* + v^) V 2 



-i + V$ 



h 



c u (h 2 ) 2 3 3 2 5(* + v / 5) V 2 
respectively. We then obtain the following result correct up to order h % : 

2i Ai 



Me u (0,h) = - 



+ 



+ 



1 + V5 3(i + y/E) V 2 
11(83^ + 89^) fh 



h\ + 65^ + 11^ fh 



2 3 3 1 5(i + y/E) V 2 
+ 0(h 7 ) 



(90) 



22335(19^-5^) V 2 . 

We note here that for general values of v not equal to an integer one obtains (up to and 
including contributions of 0{h & )) 
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Me u (0,h) = 1 + 



h 



+ 



z/ 2 + 2 



h 



+ 



v 2 - 1 \2J (z/ 2 - l)(z/ 2 -4) V2 
2(z/ 6 + 4z/ 4 - 39z/ 2 - 62) //i 



+ 0(/i 7 ) 



(91) 



(z/ 2 - l) 3 (z/ 2 -4)(z/ 2 -9) \2 t 

Next we evaluate M^(0, h) with the help of eq.(87) choosing r = and then as a check 
r = 2. In the first case we obtain the expansion 



M<p(0,h) = Mh)J u {h) - AijMJ^h) - -^-J^J^h) 

c {n ) c {n ) 

+ c MlMh)j u+2 (h) + c ^±j„ 2 (h)j v „ 2 (h) + ■■■ 



c» (h 2 ) 

In lowest orders of h 2 this is 



c» (h 2 ) 



M«(0,/i) = J (h)J v (h) + -^-±J- 2 (h)J„_ 2 (h) 

c {n ) 



which when evaluated in lowest orders of h 2 implies 



(h/2f 



.1 + 



-3 {h/2f 
2 + iV5' 2!0! 



.1 



It follows that 



Mi 1 \0,h) = -(h/2y 



-1 +i v / 5' 
2 + iV5 . 



Here, of course, v has as before the S-wave value, i.e. 



v = 2- 



(92) 



(93) 



(94) 



If we set r = 2 in eq.(87) we obtain 



— J- 2 J u+2 -J-±J, 



i/+3 



C 2 C C 4 

—^-JsJu+i H — -JqJu+a H —J-iJv + 



In lowest orders of /i 2 this is 



(95) 
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which when evaluated in lowest orders gives 



3.2 7 (fr/2)"+ 4 
1? v\ 



which reduces to 



2(z/+ l)(z/ + 2) (2 + ^5)4.3.2 

i(k/2)' r - 1 + iVr 



2 + iy/E 

in agreement with the previous result, i.e. eq.(93). 

Next we come to M^(0, h). Again we use first the method with r = in eq. (87). We 
have 



M™(0,h) = J Q (h)J- u (h) - ^-J-J^J^h) - -^-J-J^J^h) 



(96) 



In lowest orders of h 2 this is 



M { V(0,h) = J (h)J-„(h) 



c?(h* 



-J 2 {h)J- v +2{h) 



c u (h 2 ) 

which when evaluated in lowest orders of h 2 implies (with the help of the power expansion 
of the Bessel function J M (2/i)) 

\ (-i/)! +U{ ) + 2\(-v + 2)\) + 2c» Ah/2) 



and so 



It follows that 



M«(0,/ l ) = i(/ l /2) 2 (l-iv / 5) 
If we set r = 2 in eq. (87) and evaluate Afi^O, /i), we obtain 



(97) 



M^(0,/i) 



J-2J-V+2 —nJ-iJ-v+3 



—J-zJ-v+i H — z^JqJ-v+a H — —J-aJ-v + 



(98) 



27 



In lowest orders of h 2 this is 



J-2J-U+2 H — —JqJ-v+a + 



Evaluating this as before we obtain in leading orders 



Hence in leading order 



2 + iV5 



\ h / 2 f.l-( -J— 
2 l ' ' \2 + iVh 



-1/+4 



z/ + 4)! 



(99) 



M^(0,h) = (h/2f 



6 



(100) 



which is seen to be in agreement with eq. (97). With Mathematica we obtain higher order 
terms, i.e. 



\ + iy/l(h\ 2 1 + iVEfh 



6 {2) + 9 ^2 
+ ("290 + 151iVE + 120 (5 - iVt) (c 



, h 
+ ln- 

2 



+ ^ (-514 - 73iVE + 120 (5 - iVtj (c + In ^ 



(101) 



9 



+ 



^ (-290 - 151iV5 + 120 (5 + iVtj (c + In ^ 
^ (-514 + 73iVE + 120 (5 + iVtj (c + In 



(102) 



With these results we can evaluate a u and ol- v . Thus again in the dominant approximation 

1 



a u = Me u (0,h)/Ml 1 \0,h) = c»/ 



^{h/2f 



(103) 



and 



a- u = Me_„(0, h)/M { y u (0, h) = -c, 



2 + 



:(h/2f 



(104) 



It follows that 
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R M^(0,h) 
M^(0,h) 



i + y/E , -49 + 80 (C + In §) 



i + 128^(2^ + >/5) 



(-49 + 80 (C + In I)) (l5H-58 v / 5 + 120(-i + v / 5) (C + h*)) 8 
+ 737280(-7i + y/E) 

2 + iy/E (5 - 2i>/5) (-49 + 80 (C + In f)) fh 



3 360 \2 

+ ^^| 49 ( 449 + 85^) - 80 (^743 + 232iVf) (c + In | 

+ 19200 (2 + ^) (C + In ^) JQ^ (105) 

with eg = Cq 1 ' = 1 (cf.MS, p. 122, eq. (39)). 

We can now identify our quantities with those of ref. [1]. Comparison with the Dougall 
coefficients evaluated previously implies 

0^ = 1/0(^/2), a_„ = l/0(-i//2) (106) 

and so in the notation of ref. [1] 

R = <f>(-v/2)/<j>(is/2) (107) 

A remarkable feature of the expression (105) is its unit modulus, as was also observed in ref. 
[1]. It means that R is a pure phase factor 

R = e iwl ,RR* = l (108) 

When and why this behaviour occurs is discussed at the end of Appendix A. 



VIII. THE QUANTITY R IN THE GENERAL CASE 

In the general case, or for s = 1 + 2 sufficiently large so that no problems with singularities 
arise, we can evaluate R and so m£](0, h) and Me±„(0, h) by simply using the power series 
expansions of v and the standard Mathieu coefficients and, of course, the power series 
expansion of Bessel functions J^{2h). One then obtains 
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m£(o,/i) 



i fh 



[±v)\\2 



±v 



v — 1 V 2 



h 



2(z/ 2 T 3z/-7) fh 
(z/± 1) 2 (vt 1)(^ 2 -4) \2 



4(z/ 4 ^llz/ 3 -2z/ 2 ±59z/-23) ^/f 



This implies for i? 



(z/±l) 2 (z/^l) 3 (z/ 2 -4)(z/ 2 -9) V2 



1 , 2_ / h \ , 2(^ 2 +3t^7) 

1 !y 2 -l \ 2 ] (V-l) 2 (j/+l)(i/ 2 -4) I 2 



(109) 



z/!(z/ — 1)! sirniu 



2u 



1 + 



2 / \ 4 

2 I h\ 2(u 2 -3u-7) / ft 

(i/+l) 2 (iy-l)(^ 2 ~4) 1 2 



iv 2 -l \ 2 



7T 



2i/ 



1 + 



4z/ 



h 



2z/(4z/ 5 + 15z/ 4 - 32z/ 3 - 12z/ 2 + 64i/ - 111) /V 
+ (z/ 2 - l) 4 (z/ 2 -4) 2 V 2 1 



(110) 



The first few terms of the expansion of the function Me u (0, h) which is needed for comparison 
with the results of ref. [1] have been obtained in the previous section. 
From eq.(llO) we extract for later reference 



sin his 



R 



n 2 ^ 



v\(v-l)\ 



(in) 



(z/ 2 -l) 2 I 2 

This expansion will be used below in the low order approximation of the absorptivity for 
higher partial waves. 



IX. CALCULATION OF THE ABSORPTIVITY 



We consider the absorptivity in a general case, and hence allow for complex Floquet 
exponents u, which we set 

v = n + i(a + i(3) = (n - (3) + ia (112) 

where n — 2, 3, 4, ... and a and (5 are real and of 0{h A ). In evaluating the ^-matrix for small 
h A one has to be careful to make the expansions in the appropriate way. Thus we write SS* 



SS* = 



R* 



e i7 ™ \ 



(113) 
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which can be rewritten 
SS* = 



e 2 ™(l-^)(l 



R* 



1 - | cos2tt / 3.(-^ + ^)-e 2 ™ + isin27ip.(^ - 



R* 



• e WW 7 



(114) 



Here we set 



3 2i7T/3 



= l + if, e ^ a = l+g 



(115) 



where / is complex and g is real (in the S'-wave case g = — 2 -^ 2L (/i/2) 4 + 0(h 8 )). Then 



COS27T/3 = 1 - 9f/,9f/ w i(27r/?) 2 « i(sin27r/3) 2 



and 



(116) 



sin27r/5 = 



(117) 



Then 



SS* = (l + g) 



(118) 



(1-*)(1-^) 
We now consider two limiting cases. 

(i) a — > implying g — > 0. 

In this case R = R* and so 1/i? 2 ~ 0(/i 4 ). This is the case of real Floquet exponents and 

so excludes the case of S- waves. Here 



SS* 



~ 1 - 



2g/ 



l - 



^_ | sin tt/3 
R 



(119) 



since 



3/ ~ -(sin2vr/5) 2 ~ 2sin 2 7r/5. 



The absorptivity A is therefore given by 



A—l — SS* 



^_ | sin 7I-/3 



, / sin n 3 
4 



(120) 
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With the help of eq.(lll) this can be written 

Av 

4i/ 



4tt 2 (| 



A 



v\(v-l)\ 



1 + 



(u 2 -l) 2 \2 



+ 0(h 8 



47T 2 (| 



4(Z+2) 



1 + 



4is 



> 2 ~ If 



h 



0(h 8 ) 



121) 



(Z + l)!(Z + 2)! 

in agreement with ref. [1]. This can be easily evaluated (e.g. with Maple), e.g. already in 
the case of P-waves (i.e. in spite of singularities in higher order terms here omitted) and 
yields in this case 

n 2 h 12 



A 



3 2 2 



14 



(122) 



where 7 is the Euler constant (also written C). This result agrees with the result in ref. [1]. 
We observe, in particular, that logarithmic energy contributions arise in the expansion. The 
formula (121), of course, does not apply in the case of S-waves. 

(ii) (3 — > implying / — > 0. 
In this case RR* = 1, so that \R\ ~ O(h ), and we cannot expand as in the previous case. 
However, g pa 0(h 4 ), so that we can expand in powers of g. Thus 

SS*= ^ 

,1 



9 



- 2 +0(g s ) 



2-R 2 -R 
= l + 0(h 8 ) 



This is the case of complex Floquet exponents as in the S'-wave case. In this case 

SS*=l-± + 0(<?) 



(123) 



(124) 



and so 



A = % + = n2{h/2)8 + IT ( 7 ~ 12(7 + ln l } ) {h/2)U + 0(/ll6) (125) 
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in agreement with ref. [1] and the rough lowest order calculation of ref. [16]. We can also 
compute for this particularly interesting S'-wave case the amplitudes of the reflected, trans- 
mitted and incident waves A r , A t , Ai defined in Appendix A. One finds 

A r = R-± = -*^ + 0(h'), 

A t = 2% sin z/7T = — tt( h/2) 4 + 0(h 8 ) 
3 

Ai = Re- - ^- = - ^(V2) 4 + Om (126) 

We observe that in the limit h 4 — > we have A t = and A r = Ai, i.e. there is only 
reflection of the fluctuation or disturbance around the D-brane like reflection from a wall 
and no transmission, which can be interpreted as a vanishing of the disturbance on the 
brane (implying a Dirichlet boundary condition). On the basis of the analogy with the case 
of the open fundamental string between brane and antibrane in Born-Infeld theory we can 
expect that as the energy increases, transmission (i.e. absorption) sets in and becomes the 
dominant effect at high energies. This is similar to what one finds in quantum mechanics of a 
potential well of depth —V Q [23]. There the properly normalised transmission and reflection 
coefficients T(E), R(E), where E is the energy, have the behaviour T(E) — > 0,R(E) — > 1 
as E — > 0, but T(E) — ► 1 and R(E) — > as E — > oo. The high energy behaviour of the 
effect considered here can presumably be investigated with the help of laxge-h asymptotic 
expansions of modified Mathieu functions which we expect to be formally (i.e. apart from 
sign and complex i changes) similar to those of periodic Mathieu functions with a parameter 
q defined as the solution of 

(Z + 2) 2 = -2h 2 + 2hq + 0(h°). (127) 

The Floquet exponent is then given by [12] (p. 210), [25] 

- Ah r l + 3(g 2 + l)/64/i ^/log/i 



Tie 

COS 1XV + 1 



(8h)i/ 2 



r((l-g)/4)r((3-g)/4) + °V ^ 



(128) 



One observes that again logarithmic contributions in the energy appear. 



33 



X. CONCLUSIONS 



In the above we considered the impingement of a massless scalar field on a D3 brane 
in 10 dimensions and calculated the S matrix and partial wave absorption and reflection 
amplitudes and rates for this process. Instead of coefficients introduced by Dougall for the 
expansion of the modified Mathieu functions involved, we used (in the low energy domain) 
rapidly convergent series in terms of products of Bessel functions. We demonstrated that the 
Mathieu function coefficients are such that many different expansions in terms of products of 
Bessel functions all yield the same low energy power series for the modified Mathieu functions 
of the first kind. We think, this is the best way to evaluate the absorption rates of the problem 
in the low energy domain. The leading term matching procedures of refs. [15,16] maybe select 
dominant terms of the expansions considered here. Since the metric considered is extremal, 
one can visualise the absorption of the partial waves of the scalar field as absorption into 
the brane or black hole with vanishing event horizon (examples with nonvanishing horizon 
have for instance been treated in [24]). Since several other string models lead also to the 
modified Mathieu equation in analogous contexts, the above considerations, which have 
definite advantages over those involving the coefficients of Dougall, may be of wider interest. 
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Appendix A 

Here we recapitulate the main steps of the derivation of the S-matrix. We follow ref . [8] , 
but instead of repeating the steps there, we emphasize those which have not been written 
out explicitly there. For ease of comparison we consider the repulsive potential which means 
simply that the considerations below employ coupling g as in ref. [8] instead of g used above. 
The two cases are trivially related through 

go = ig (A.i) 

In the repulsive case we have a regular solution y reg of eq.(5)at r = 0, i.e. one proportional to 
exp(—g I r) . The variable of the cylindrical functions involved in (z, h) is uo = 2h cosh z = 
(ig/r + kr). Thus in leading order for small h 2 and r — > 0(z — > — oo) we can write 

y reg = r^M?Xz,h) ^~-°V/2 

(A.2) 

If we let Rez — > — oo here and then replace z by— z, the solution has the asymptotic be- 
haviour e tkr . The series expansion defining the Bessel function J u has the following important 
property for integers n 

J u (2hcosh(z + inn)) = exp(invir)J u (2hcoshz) (A. 3) 

so that 

M^\z + inn, h) = exp(invn)M^ (z, h). (A.4) 

Since for Me u (z, h) also 

Me v {z + inn, h) = exp{invn)Me v {z, h) (A. 5) 
we have the proportionality 

Me v (z, h) = a v (h)Ml l \z, h) (A.6) 



HW(u) + 0{h 2 ) 



ng 



e r e 



l + 0(h 2 ) 
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with (e.g.) 



a v {h) 



Me„(0, h) 



(A.7) 



As mentioned earlier, expansions in terms of cylindrical functions like (85) converge uni- 
formly only in domains |cosli2;| > 1, whereas the expansion (32) of Me v (z,h) converges 
for all finite complex values of z. Hence we match M^\z,h) in the domain $lz < to a 
linear combination of M„ 3 \z, h) and Mu 4 \z, h) in the domain $lz > by matching both to 
a combination of Me u (z, h) and Me_„(z, h) in the intermediate domain. We have 



z = log \fkj~gr (+) in/A. 



(A.8) 



In the domain of r close to zero we set 



d 
dr 



r V2 M G») = r i/2 ^Me^ + (SMe^j 



r V2 M (3) 



d 

a— I r J 

dr 



(A.9) 



where a and (3 have to be determined. In the domain of large r we set, with constants 
a', (3', A, B, which have to be determined 

r 1 ' 2 (a'Me v + = r 1 ' 2 (aM® + BMf'j , 

4r (^) + 4 t V2 ^) ^ 4 (^) + *3F (^) • (A-10) 

We match the Me„, Me_„ combination (variable z) on the left to that on the right (variable 
— z) at sftz = 0(r = y/g/h), so that 



r 1/2 (aMe u + /3Me_J = r 1/2 (a'Me v + /3'Me_„ J 



d 

a— | r 

dr 



l/2 Me)j +pj-(rV 2 Me_ 



z=+iw/4 



z=—in/4 



(All) 



z=— J7r/4 



Since Me^z) = Me_„(-z) and at r = ^^/gjk,z = ±iir/A, also 



dr 



k\ 1/2 d 



g) dz' 
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the latter become 



&Me_ v + (3Me v = a'Me u + /3'Me_„ 

a^-Me_ v + p-^-Me u = d^-Me v + f3'-^Me_ v (A.12) 
dz dz dz dz 

From these equations we obtain immediately 

d = P, (3' = a (A. 13) 

From eqs.(A.9) we obtain (W meaning Wronskian) 

= W\Mf\Me. v \ = W[M<?\Me u ] 

W[Me v ,Me- v ]' W[Me v ,Me- v Y 1 ' ' 

From (A. 10) we obtain similarly 

-W[M„ 3 \ Me u }W[Me u , A/£ 4) ] + W[M^ , Me„ u ]W[Me_ u , il#] 



A = 



W[M^\Mi 4) ]W[Me v , Me.,] 



B = W[Mj?\Me u ]W[Me v , M< 3) ] - W[M^ , Me. u }W[Me. u , M^ 3) ] 

WlMl 3) ,M^}W[Me u ,Me_ u ] 

We now use eq. (A.6) and Wronskians W[M^\m^\ = given in ref. [12], i.e. 

[3,4] = --, [1,3] = -[1,4] = - (A.16) 

71 71 

and the circuit relation ( [12], p. 169) 

M (i) = e ivn M {i) _ i sin U7r M^ . (A. 17) 

Then 

W[Me u , Me- V ] = _ 2sm u% av0L _^ 

71 

W[Me. u ,M^} = -e~ i,J *a- v , W[Me- v ,M^\ = --e il/7T a. u , (A.18) 

7T 7T 

With these expressions A and B are found to be 

f ( OLy OL — u \ t-, 1 ( Otv — 2iv7Y 



2ismv7i\a- u a v J 2i sin utt \ a- v a v 



5 = — - JZ__ e -^r=!l (A.19) 
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The regular solution thus continued to r = oo is then 

AM^(z,h)+BM^(z, h) 
2 \ 1/2 f 



t.1/2 
yreg — ' 



(JL^j | Ae ikr e- i{l/+1 ^ +e-^Be- ikr e^ + ^ (A.20) 



In terms of the variable z and with R = a u /a_ u this can be written 



2r 



2/i7r cosh 2 



1 



2h% cosh z / { R 

-tun 

+ i(Re tun - l_) e -2^cosh* l ^ A 21 ^ 

If we take A^ = (Re lU7T — 6 ^ ) as the amplitude of the incident wave, the amplitudes A r and 
A t of the reflected and transmitted waves are A r = R — and A t = 2i sin z/7r respectively in 
agreement with ref . [1] . With the definition of the S'-matrix in the partial wave expansion of 
the scattering amplitude f(9), with x = cos 9, for n space dimensions (here we have n — 6) 
given by [21] 

ikr i °° (- ^ 

e ^ + W)^=w - 2( ^ r) (n-i)/ 2 E + ("l)'*"- 1 ^ j^(c™ fl) (A.22) 



where 



VticosO) = \f^T/ 2 - l T(n/2 - + J - 1)C,( 

V 7T Z 



cost 



and C;(cos^) is a Gegenbauer polynomial, we obtain for this 

R — — 

S = «_. e -*T' (A.23) 

It is easy to verify that for v real and R = e y real, unitarity is preserved, i.e. unity minus 
reflection probability = transmission probability, i.e. 

I-R- |2sin utt\ 2 

We observe that this relation remains valid if the real quantity R = e y and the pure phase 
factor e mv exchange their roles, i.e. if R becomes a pure phase factor and e mv a real 
exponential. The latter is precisely what happens in the S-wave case of the attractive 
potential discussed above. 
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Appendix B 

Below we give the explicit form of the first three terms of the small-ft 2 perturbation 
expansions of the leading coefficients c 2 " of expansions of modified Mathieu functions in the 
S-wave case (1 = 0). The nonleading terms have been obtained with Mathematica. 



c 2 v _ i + V$ f h Y + U(5 -W5) fh\ 6 + (7231H- 17746^5) AV^^iij 



Co" + 2 3 5(-2* + v^) 2 V 2 / 2 7 3 3 5(-2^ + v^) 3 V 2 



cZ u 2 _ 1 (hy 3 + lGiy/E fhy + 5146 - 391i^5 fh\ , /1 , ; ,u, 



3 V 2 / 23 3 3 V 2 / 26355 ~ V 2 



cl u _ i - V5 + ll(5i + VE) fh\ 4 + 29QQV5 - 28889i / h\ , , ; 



c " i + 2 3 3 1 5(i + >/5) V 2 / 2 7 3 5 5(2i + V 2 



CI4 J_fh) 18 + 125^^5 fhy 1303^- 25774 f h\ 13 

77 2^3 v^y + 26335 V 2 / + 2^ V 2 ' 



,n "A 6 39559i- 4597^5 //A 10 m 



cT" 1 fhy 290^ + 49^ (h 



c Q v 2 + i v / 5V 2 / 2 2 3 3 5(i + y/E) V 2 



dg 1 //A 6 19 + 157ty^5 //A 10 132035 - 5159^ / /A 14 15 

^ ~ ~23p5 V2 J 2^52 V 27 + 2 7 3 6 5 3 7 I2J + ( 
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